Dynamic instabilities and memory effects in vortex matter 
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Understanding the nature of flow is essential for the resolution of a wide 
class of phenomena in condensed matter physics, ranging from dynamic friction, 
through pattern formation in sand dunes, to the pinning of charge density waves. 
The flux line lattice in type II superconductors serves as a unique model system 
with tunable dynamic properties. Indeed, recent studies have shown a number 
of puzzling phenomena including: (i) low frequency nois (ii) slow voltage 
oscillations!'!, (m) history dependent dynamic responsdzHll, (iv) memory of the 
direction, amplitude, duration, and even the frequency of the previously applied 
currentlll'B, (v) high vortex mobility for ac current with no apparent vortex 
motion at dc currentllUll'Ill, and (vi) strong suppression of an ac response by small 
dc bias§. Taken together, these phenomena are incompatible with the current 
understanding of vortex dynamics. By investigating the current distribution 
across single crystals of 2H-NbSe2 we reveal a generic mechanism that accounts 
for these observations in terms of a competition between the injection of a 
disordered vortex phase at the sample edges, and the dynamic annealing of this 
metastable disorder by the transport current. For an ac current, only narrow 
regions near the edges are in the disordered phase, while for dc bias, most of the 
sample is filled by the pinned disorder, preventing vortex motion. The resulting 
spatial profile of disorder acts as an active memory of the previous history. 

In conventional superconductors like NbSe 2 the anomalous phenomena are 
in the vicinity of the 'peak effect' where the critical current I c increases 
sharply below the upper critical field H c2 , as described in Fig. la. The peak effect marks 
a structural transformation@'!IHi30 of the vortex lattice: Below the peak region an ordered 
phase (OP) is present, which is dominated by the elastic energy of the lattice and is, there- 
fore, weakly pinned. On approaching the peak region, however, the increased softening of 
the lattice causes a transition into a disordered vortex phase (DP), which accommodates 
better to the pinning landscape, resulting in a sharp increase in I c . In high-temperature 
superconductors like Bi 2 Sr 2 CaCu208, this situation is equivalent to the second magnetiza- 



tion peaklll, where the ordered Bragg-glass phase is believed to transform into a disordered 
soliofHl Figure la shows the I c measured at various frequencies. On the high temperature 
side of the peak effect I c is frequency independent; in this region the DP is thermodynami- 
cally stable. In contrast, on the low temperature side, a significant frequency dependence is 
observedlliJii; in this region all the unusual vortex response phenomena appearltfl&H. 
As described below, a dynamic coexistence of the OP and a metastable DP is established in 
this region in the presence of an applied current. We first outline the proposed mechanism, 
and then present the experimental evidence. 

The first important ingredient of the proposed model is the observation that in NbSe2 
the DP can be readily supercooled to below the peak effect by field cooling, where it remains 
metastable, since the thermal fluctuations are negligiblelHS. This supercooled DP is pinned 
more strongly and displays a significantly larger critical current J^ s as compared to J° rd 
of the stable OP. An externally applied current in excess of J dls serves as an effective tem- 
perature and 'anneals' the metastable DP as observed by transporti, magnetic response^, 
decorationil, and SANS experiments^ on NbSe2- 

The second ingredient of the model is the presence of substantial surface barriers@, as 
observed recently^ in NbSe2. Consider a steady state flow of an OP in the presence of 
a transport current. In the standard platelet strip geometry, in a perpendicular field B, 
vortices penetrate from one edge of the sample and exit at the opposite edge. In the absence 
of a surface barrier, vortex penetration does not require any extra force. As a result, the 
vortices penetrate close to their proper vortex lattice locations, as dictated by the elastic 
forces of the lattice. In the presence of a surface barrier, however, a large force is required for 
vortex penetration and exit, and hence much of the applied current flows at the edges in order 
to provide the necessary driving forceii~il. The surface barrier is known to be very sensitive 
to surface imperfections. Therefore, the penetrating vortices are injected predominantly at 
the weakest points of the barrier, thus destroying the local order and forming a metastable 
DP near the edge, which drifts into the sample with the flow of the entire lattice. (Note 
that steps on the surface or extended defects in inhomogeneous samples could also act as 



injection points of the DP). The applied current, therefore, has two effects: the current that 
flows at the edges causes 'contamination' by injecting a DP, while the current that flows in 
the bulk acts as an annealing mechanism. The observed dynamic instabilities and memory 
phenomena arise from the fine balance between these two competing processes. 

The annealing process is sensitive to the exact location on the H — T phase diagram. 
Below the peak effect, the DP is highly unstable and therefore its relaxation time r r in the 
presence of a driving force is very short. As a result, it anneals rapidly over a characteristic 
'healing' length L r = vr r , where v is the vortex lattice drift velocity. The corresponding 
profile of the local critical current J c {x) should therefore decay from J d%s to J° rd over the 
characteristic length scale L r , as illustrated by the dotted line in Fig. lb. Note that L r 
and r r are generally current dependent and decrease dramatically at elevated currents!!. On 
the other hand, near the peak effect the free energies of the DP and OP are comparable 
and therefore the 'life time' of the disordered phase r r and the corresponding L r are very 
large. As a result, the front of the DP, given by xj(t), progressively penetrates into the 
bulk as shown by the solid line in Fig. lb, until the entire sample is contaminated. In this 
situation the experimental, steady state dc critical current I dc does not reflect an equilibrium 
property, but rather a dynamic coexistence of two phases. It is given by I dc = d J c (x)dx ~ 
dL r Jf s + d{W - L r )J° rd for L r < W, and if ~ If 3 = dWJf s for L r > W, where d and 
W are the thickness and width of the sample (neglecting, for simplicity, the surface barrier 
edge currents, and assuming, for example, an exponential decay of J c (x) in Fig. lb). 

The ac response of the system should be distinctly different since the contamination 
process occurs only near the edges, where the disordered lattice periodically penetrates and 
exits the sample. For a square wave I ac of period T ac = 1/ f, by the end of the positive half 
cycle the DP occupies the left edge to a depth of x^f, as illustrated by the solid curve in Fig. 
lb. During the negative half cycle the DP on the left exits the sample, while a DP on the 
right edge penetrates, until at t — T ac a mirror-image profile is obtained, as shown by the 
dashed line. Assuming x^ c < W < L r , the effective I c observed by ac transport measurement 
is given by 7" c ~ dx°£J dts + d(W — x c £)J° rd . Thus, an ac current necessarily contaminates 
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the sample less than a dc current of the same amplitude, and therefore < 1^ always, as 
seen in Fig. la. In addition, since x°£ decreases with frequency, J? c should decrease with / 
explaining the frequency dependence of in Fig. la. Furthermore, and most importantly, 
at sufficiently high frequency J" c should approach the true I c of the stable phase. The steep 
increase of the 881 Hz I® c data in Fig. la (open circles) therefore indicates that the OP 
transforms sharply into the DP at the peak effect. In contrast, the smooth behavior of 
reflects rather the dynamic coexistence of the two phases in which L r gradually increases 
and diverges upon approaching the peak effect from below. From Fig. la we can evaluate 
L r and r r . For example, at T=5.1K, Ifr ~ 50 mA is about half way between I° rd ~ 5 
niA and the extrapolated lf s ~ 100 mA, which means that L r ~ 0.5W 7 = 170/im. The 
was measured at a voltage criterion of 1 /iV, which translates into vortex velocity of 
i;~4x 10~ 3 m/sec, and hence r r = L r /v ~ 4 x 10~ 2 sec. This value is well within the range 
of the relaxation times measured previously! by applying a current step to the field-cooled 
matastable DP. 

We now provide a direct experimental manifestation of the key aspect of the model, 
which is the spatial variation of the disorder and J c (x), and of the transport current distri- 
bution that traces this J c (x) (see Fig. lb). We have used Hall sensor arrays to measure the 
ac transport current self-induced fieldiHl B ac (x), which is then directly inverted into the 
current density distribution J ac {x) using the Biot-Savart law, as described previously^! (see 
inset to Fig la). Figure 2 shows the corresponding current profiles J ac {x) measured at dif- 
ferent frequencies. At high /, the DP with the enhanced J c is present only in narrow regions 
near the edges (481 Hz data). As the frequency is reduced, x^f grows and correspondingly 
the enhanced J ac (x) flows in wider regions near the edges. Note that our measurement 
procedure (see Fig. 2) provides the time averaged local amplitude of J ac (x), which is much 
smoother as compared to the sharp instantaneous profiles in Fig. lb. 

We confirm the above finding independently by measuring the corresponding ac resistance 
of the sample R ac {f) as shown in Fig. 3a. At high frequencies most of the sample is in the 
low pinning OP and therefore R is large. As / is decreased, progressively wider regions 



near the edges become contaminated with the more strongly pinned DP and thus R ac {f) 
decreases. If the applied I ac is larger than I^ c , a finite R will be measured at all frequencies, 
however, if 7" c < I ac < (see Fig. la) the measured R will vanish as / — * 0, as observed 
in Fig. 3a. This explains the surprising phenomenon of finite vortex response to ac current, 
while for dc drive the vortex motion is absent0il0. From R ac {f) one can directly calculate 
the width of the disordered regions by noting that x d c equals the distance the entire lattice 
is displaced during half an ac period, x d c (f) = v/2f = R(f)I ac /2fLB, where L is the 
voltage contact separation. The open circles in Fig. 3b show x a d obtained from R ac (f), 
while the open squares show the x d c derived directly from the J ac {x) profiles of Fig. 2. 
The good correspondence between the two independent evaluations of x d c demonstrates the 
self-consistency of the model. 

Next we address the extreme sensitivity of the ac response to a small dc biasEl, as shown 
in Fig. 4a, where R ac is presented as a function of a superposed Idc- A dc bias of only 10 to 
20% of I ac suppresses R ac by orders of magnitude. This behavior is a natural consequence 
of the described mechanism since the dc bias contaminates the sample very similarly to the 
pure dc case, except that L r is now renormalized as following. For I dc <C I ac , the vortices 
move back and forth during the ac cycle, with a forward displacement being larger by about 
2Idc/I a c- Therefore, a vortex that enters through the sample edge and reaches a position x, 
accumulates a much longer total displacement path of xl ac /ld c - Since the annealing process 
of the DP depends on the total displacement regardless of the direction, the lattice at this 
location is thus annealed substantially, as if the effective L r is reduced to L e J$ ~ L r I dc /I ac . 
Thus, at very small biases, the DP is present only within x a d c from the edges, as in the 
absence of a bias, where the disordered vortices exit and re-penetrate every cycle. Vortices 
that drift deeper into the bulk under the influence of I dc are practically fully annealed due 
to the very short L e J* . As a result in Fig. 4a the initial decrease of R ac up to I dc ~ 2 
mA is relatively small. The corresponding J ac {x) in Fig. 4b at I dc = 1.7 mA shows narrow 
contaminated regions near the edges, very similar to the zero bias case in Fig. 2. However, 
as Idc is increased, L e J* grows and the bulk of the sample becomes contaminated by the 
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penetrating DP, leading to a dramatic drop of R ac . In this situation, J ac (x) at Id c = 5.7 
mA shows a wide region of DP at the left edge. When is inverted to -5.7 mA, a similar 
situation is observed, but now the vortices and hence the DP penetrate from the right edge, 
as expected. 

The revealed mechanism readily explains a wide range of additional reported phenomena: 
{%) The history of the previously applied current is encoded in the spatial profile of the lattice 
disorder, which is preserved while the current is switched off due to negligible thermal 
relaxation. Upon reapplying the current, the vortex system will display a memory of all the 
parameters of the previously applied current including its direction, duration, amplitude, 
and frequency, as observed experimentallyfilE. {%%) Application of a current step / < 
to a sample in the OP, results in a transient response which decays to zero since the DP is 
able to penetrate only a limited distance. The resulting new I c of the sample is given by 
the condition that I c = I, as derived by fast transport measurementSEl. Such transient 
phenomena, would also display characteristic times shorter or comparable to the vortex 
transit time across the sample, in agreement with observations00. (Hi) The competition 
between the contamination and annealing processes is expected to result in local instabilities 
causing the reported noise enhancement below the peak effect!! (see also Fig. 4a). (iv) 
Related phenomena should be observed in high-temperature superconductors in the vicinity 
of the peak effect associated with the melting transition, or near the second magnetization 
peak, consistent with experimentsi"i'Ei@. (v) In high-temperature superconductors there is 
an additional consideration of thermal activation of vortices over the surface barriers, which 
may explain the reported slow voltage oscillationslHl If the thermal activation rate is higher 
or comparable to the driving rate, the slowly injected lattice will be ordered in contrast to 
the DP injected at higher drives. Thus, at a given applied current, if the bulk of the sample 
is in the OP, much of the current flows on the edges, rapidly injecting a DP through the 
surface barrier. Once the bulk gets contaminated, the resulting slower vortex motion causes 
again injection of an OP. This feedback mechanism can explain the voltage oscillations!'! 
in YBa2Cu30y and similar narrow band noise!'! in YBa2Cu30y and B^S^CaC^Os with 



characteristic frequencies comparable to the inverse transit time, (vi) Finally, the described 
phenomena should be absent in the Corbino disk geometry where vortices do not cross the 
sample edges. Our studies of NbSe2 in this geometry confirm this prediction, as will be 
published elsewhere. 
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FIGURE CAPTIONS 

Fig. 1. The experimental setup (inset), the critical current vs. temperature in the 
vicinity of the peak effect (a), and (b) a schematic plot of the dynamic coexistence of the 
ordered phase (OP) with a metastable disordered phase (DP). 

Experimental. Several Fe (200 ppm) doped single crystals of 2H-NbSe 2 were investigated. 
Here we report data on crystal A of 2.6 x 0.34 x 0.05 mm 3 and T c = 6. OK, and crystal B of 
2.4 x 0.29 x 0.02 mm 3 with T c = 6.05K. Four electrical contacts were attached to the top 
surface for transport measurements, with the voltage contact separation of 0.6 ± 0.2 mm. 
The bottom surface of the crystal was attached to an array of 19 2DEG Hall sensors 10 x 10 
fim 2 each (inset). The vortex lattice was initially prepared in the OP by zero-field-cooling 
to a low temperature at which the DP is unstable, and then slowly heated to the desired 
temperature in the presence of a constant field H a applied parallel to the c axis. 

(a) The peak effect in critical current I c in NbSe 2 crystal A vs. temperature, as measured 
with a dc drive, lj c (□), and ac drive, J? c , at 181 (■) and 881 Hz (Q). The critical current 
was determined resistively using a voltage criterion of 1 /iV. At low temperatures only the 
stable OP is present. On the lower temperature side of the peak effect a metastable DP 
coexists dynamically with the stable OP resulting in a frequency dependent J c . On the 
high temperature side of the peak effect only the stable DP is present with no anomalous 
behavior. 

(b) Schematic plot of the local critical current density J c (x) across a crystal of width 
W. j£ ls and J° rd are the values of the critical current density in the fully disordered and in 
the OP, respectively. For L r <^ W the DP relaxes rapidly into OP resulting in the dotted 
J c {x) in the steady state dc flow. For L r > W the DP penetrates to a depth Xd(t) following 
the application of dc current at t — (solid curve). For an ac current at t = T ac /2 the DP 
occupies the left edge to a depth of x a d c (solid curve), and symmetrically the right edge at 
t = T ac (dashed curve). 

Fig. 2. Current density profiles J ac (x) in crystal B obtained by inversion of the self- 
induced field measured by the Hall sensors. Shown are three frequencies / = 22 (•), 181 
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(□), and 481 Hz (■). A square wave ac current I ac was applied and the corresponding self 
induced magnetic field B ac (x) across the crystal was measured by the Hall sensors using a 
lock-in amplifier (see Fig. 1 inset). By using the Biot-Savart law the B ac (x) was directly 
inverted^ into the current density profiles J ac (x). The width x d c of the highly pinned DP 
near the edges grows with decreasing frequency as expected. The measured J ac (x) is the 
magnitude of the local current density averaged over the ac cycle period. As a result, J ac (x) 
reflects a time-averaged superposition of solid and dashed J c (x,t) profiles in Fig. lb, which 
are present separately during the positive and negative half-cycles. Close to the edges the 
high J is present most of the time, while close to x d c it is present only a small fraction of the 
ac period as the DP moves in and out of the sample. Therefore, the time-averaged J ac (%), 
decreases smoothly from the edge to x a d . A more detailed analysis shows that the first- 
harmonic measurement by the lock-in amplifier Fourier transforms the sharp instantaneous 
J c (x,t) of Fig. lb into the observed smooth J 6 ac (x) = J° rd + {Jf s - J° rd )(l + cos(-Kx / x a d c )) /2 + 
(Jf s — J° rd )(l-cos(n(x + x d c -W)/x d c ))/2. The second and third terms hold at < x < x a d c 
and W — x a d c < x < W, respectively, and are zero otherwise. 

Fig. 3. Frequency dependence of the resistance R ac (f) and of the width of the disordered 
regions x a d in crystal B. (a) At high frequencies R ac is large since most of the sample is in 
the weakly pinned OP. As the frequency is decreased the disordered regions increase and R ac 
drops sharply, when x d c reaches values close to the sample width. In the limit of zero / the 
resistance is zero since the applied current is lower than I dc . (b) The corresponding width 
of the DP near the edges x d c derived from the R ac {f) data (O) an d from the J ac {x) current 
profiles of Fig. 2 (□). As expected, x a d increases monotonically with decreasing frequency. 
The drop of the x a d data (O) & t very low frequencies is an artifact. At such frequencies the 
instantaneous vortex motion is present mainly at the onset of the square wave I ac pulses, and 
decays towards zero during the pulsei. In this situation the first-harmonic R ac measurement 
by the lock- in amplifier underestimates the integrated vortex displacement. In addition to 
the frequency dependence, x d c also changes significantly by varying the amplitude of the ac 
current I ac . Here 7 ac =10 mA was chosen such that x a d becomes comparable to the sample 
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width W at low frequencies. 

Fig. 4. Measured resistance R ac at I ac = 20 mA as a function of dc bias I dc (a) and the 
corresponding distribution (b) of the ac current J ac (x) in crystal A. (a) At low bias Idc < 
2 mA only the regions near the edges are contaminated by the DP since L e Jf is very short, 
resulting in only a moderate decrease of R ac . For I dc > 2 mA the contamination becomes 
substantial and the significant decrease of R ac is accompanied by enhanced noise due to 
the local instabilities during the competing contamination and annealing processes. At still 
larger dc bias most of the sample becomes contaminated and R ac drops below our noise 
level. The corresponding J ac (x) profiles in (b) show that for positive Idc = +5.7 mA (•) a 
substantial part of the sample is contaminated from the left edge where the vortices enter 
into the crystal, and similarly from the right edge for negative bias I dc = —5.7 mA (■). 
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